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ABSTRACT 



Aims. The presence of depletion (freeze-out) of CO around low-mass protostars is well established. Here we observe 84 
candidate young high-mass sources in the rare isotopologues C 17 and C 18 to investigate whether there is evidence for 
depletion towards these objects. 

Methods. Observations of the J = 2 — > 1 transitions of C ls O and C 17 are used to derive the column densities of gas towards 
the sources and these are compared with those derived from submillimetre continuum observations. The derived fractional 
abundance suggests that the CO species show a range of degrees of depletion towards the objects. We then use the radiative 
transfer code RATRAN to model a selection of the sources to confirm that the spread of abundances is not a result of 
assumptions made when calculating the column densities. 

Results. We find a range of abundances of C 17 that cannot be accounted for by global variations in either the temperature 
or dust properties and so must reflect source to source variations. The most likely explanation is that different sources show 
different degrees of depletion of the CO. Comparison of the C 17 linewidths of our sources with those of CS presented by other 
authors reveal a division of the sources into two groups. Sources with a CS linewidth > 3 kms -1 have low abundances of C 17 
while sources with narrower CS lines have typically higher C 17 abundances. We suggest that this represents an evolutionary 
trend. Depletion towards these objects shows that the gas remains cold and dense for long enough for the trace species to 
deplete. The range of depletion measured suggests that these objects have lifetimes of 2 — 4 x 10 5 years. 
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1. Introduction 

Many questions regarding the formation mechanism of 
high-mass stars remain unanswered. Young high-mass 
stars form within massive cores and display features com- 
mon to low-mass star formation such as outflows and 



in so me cases disks (e.g. IShepherdl l2005t ICesaroni et al 



20071) . However the different stages through which a core 
forming a high-mass star evolves remains unclear. 

In the cold, dense conditions within cores prior to their 
collapse, trace molecules freeze-out onto the dust grains 
forming icy mantles. This process reduces the abundances 
of species in the gas phase until the grains are heated, 
desorbing the molecules off the grains and returning them 
back into the gas. This freeze-out and evaporation cycle 
is potentially an important indicator of the age of the 
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material in star forming cores and how the material has 
evolved. 



The transitions of the two rare isotopologues of CO, 
C 18 and C 17 which have abundances with respect 
to H 2 of ~ 2 x 10~ 7 and ~ 5 x 10~ 8 respectively 



(jFrerking. Langer fc Wilsonlll982D , ar e often employed to 
probe the inner reg i ons of dense cor es jTafalla et alJl2002t 



Casclli etaJJ 1199 



3giom 
9; R. 



cdmanetal 



2002ri 



as their tran- 
sitions are typically optically thin, especially those of 
C 17 0. A number of studies of various, but mainly low- 
mass, star forming regions using these species have pro- 
vided evidence for depletion of CO towards th e cen 
tres of cores (Caselli et al.l Il999t Kramer et al.l 1999; 



Tafalla et al1l2002HWillacv et al.lll998HSavva et al JI2003 ) 
with the degree of depletion varying from from factors 
of a few to over an order of magnitude compared to the 
canonical CO abundance. 
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In this paper we present new observations of candidate 
high-mass star forming regions in C 17 and C 18 made 
in order to look for evidence of depletion in these regions. 



2. Observations of C 17 and C 18 



•18| 



Sridharan etal] (|2002l hereafter SBSMW) searched the 
IRAS catalogue and identified 69 point sources which rep- 
resent potentially massive, deeply embedded protostars 
in the Galactic plane. These s ources have been studied 
extensively by iBeuther et all (|2002al lbT) and have been 
mapped with SCUBA, the James Clerk Maxwell Telescope 
(JCMT) bolometer array, a t 850 ^m and 450 /im by 
Williams. Fuller fc Sridharan! (|2004l . hereafte r WFS04) as 
well a s being surveyed for evidence of infall bv I Fuller et al.l 
( 2005h . In many cases the submillimetre and millimetre 
continuum maps show more than a single peak indicat- 
ing multiple sites of potential star formation. In total 112 
850 /im peaks were identified of which we have observed 
84 in C 17 0. 

We observed the J = 2 — > 1 rotational transitions of 
C 18 (219.560 GHz) and C 17 (224.714 GHz) towards 
31 sources during May 2004 at the JCMT0 in Hawaii as 
part of observing program M04AU47. We used the RxA3 
receiver with the DAS autocorrelator with channels of 
78 kHz in position switching mode using a reference po- 
sition of (+600", +600"). The reference positions were 
chosen to be free of C 18 emission having been checked 
and used for previous molecular observations. The point- 
ing was checked every couple of hours. The typical zenith 
opacity at 225 GHz was between 0.15 and 0.4 and we 
achieved a typical rms of <0.25K. The data were base- 
line subtracted using the SPECX package following nor- 
mal routines. 

As we are primarily interested in deriving column den- 
sities, in further observations we chose to concentrate on 
the C 17 line which has lower optical depth. A further 
53 sources were observed in C 17 J = 2 — > 1 in July and 
August 2005 (in observing program M05BU46). In ad- 
dition we observed a selection of these sources in the 
J = 3 -> 2 transition of C 17 (337.061 GHz) with the 
RxB3 receiver. This time we achieved a typical rms of 
<0.15K. These data were processed in the same way as 
the earlier data. 

As the telescope beam is not perfectly coupled to the 
source it is necessary to apply a correction factor to con- 
vert antenna temperature into a main beam temperature 
where 

T — A 
TJmb 

At the JCMT the main beam efficiencies are 0.69 at 
220 GHz with a beam size of 20" and 0.63 at 337 GHz with 



1 The James Clerk Maxwell Telescope is operated by The 
Joint Astronomy Centre on behalf of the Particle Physics 
and Astronomy Research Council of the United Kingdom, 
the Netherlands Organisation for Scientific Research, and the 
National Research Council of Canada 



a beam size of 13". The positions and dates of observations 
for all the sources are given in Table [5] 

3. Analysis 

3.1. Line Profiles 

The line parameters were initially measured using SPECX 
by fitting Gaussian line-profiles to derive the source veloc- 
ity (vi sr ), line width (Av), peak line flux (T\) and inte- 
grated line intensities. 

All the C 18 data were fitted with a single Gaussian. 
These fits produced typical peak temperatures of T\ ~ 
2.5 — 7.5 K, and linewidths of Av ~ 1.5 — 3.5kms _1 . 
However examination reveals that a significant fraction 
of these C 18 data (~65%) are distinctly better fit, with 
smaller residuals, by the sum of two different velocity com- 
ponents. The results of two components fits are also given 
in Table [7] Those sources with multiple components can 
be divided into two categories; those where the 2 com- 
ponents are of approximately equal width, but are offset 
in velocity, and those where there is a definite separation 
into a broad and a narrow component, with the broad 
component possibly being a related to the outflow from 
the source. In this latter case the typical linewidths of the 
broad and narrow components are Auf, r > 3kms _1 and 
Av n < 2kms~ 1 . 
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Fig. 1. Example spectra of C 18 line profiles for (a) WFS30 
and (b) WFS29 which are best fit by 2 Gaussian components 
shown as the smooth curve over each line. 



The integrated intensities for C 17 have been calcu- 
lated from a Gaussian fit using SPECX. However unlike 
the C 18 0, the transitions of C 17 have hyperfine struc- 
ture. Simple Gaussian fits to the C 17 data therefore 
overestimate the intrinsic velocity dispersion. The C 17 
line parameters were therefore fitted with the known hy- 
perfine structure (hfs) of the respective transitions us- 
ing METHOD HFS in the CLASS package assuming all 
components have equal excitation temperatures. For the 
J = 2 — » 1 transition there are 9 hyperfine components, 
while there are 14 for the J = 3 — ► 2 transition. The line 
parameters, the peak intensity, the vi sr , the lincwidth and 
the integrated intensity are given in Tables Q] and [2] The 
optical depths calculated from the hfs fits imply these 
sources are optically thin in C 17 with r ~ 0.1. 
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Table 1. Observed line parameters. The columns give the peak intensity, velocity, linewidth and integrated intensity for C 18 
and C ir O respectively. For C 17 the first three parameters are derived from hyperfine fitting of the data using CLASS. The 
integrated intensity was obtained using a single Gaussian fit. 



3.2. Ratio of Integrated Intensities 

The ratio of C 18 to C 17 line intensities can constrain 
the opt ical depth of the line if the abundance ratio is 
known ( Ladd. Fuller fc Deanel 1998 ). Penzias conducted 
a survey of 1 5 massive star forming regions lying in the 
galactic disk ( Penziasl 1981 ) and derived a value of R = 
[ 18 0]/[ 17 0] = 3.65 ± 0.15 and found no gradient based 
on galactocentric distance. Further work by different au- 
thors have produced values, assuming that [ C 18 Q1/[C 17 01 
= [ 18 0]/[ 17 0], rangin g from R = 2.9 ± 1.2 (ISheffer et al 



20021) to 4.1 5 ± 0.52 l|Bensch et al.ll200lh . More recently 
Laddl ( 20041 ) conducted a survey of 648 lines of sight to- 
wards 5 star forming regions in the Taurus molecular cloud 
and combining the results for the individual regions con- 
cluded that R = 4.0 ± 0. 5, a value consistent with an 
analys is of other clouds by IWouterloot. Brand &: Henkel 
( 20051) . 

A comparison of our measured C 17 and C 18 inte- 
grated intensities is shown in Figure [2] Assuming equal 
excitation temperatures and beam filling factors for the 
two species, if the emission is optically thin one would ex- 
pect the ratio of the integrated line intensity of the two 
species to equal the abundance ratio. A standard value 
of 3.65 for the abundance ratio of C 18 to C 17 is indi- 
cated by the solid line. Of the sources, 5 lie approximately 
on the line implying that for the C 18 and C 17 lines 
are optically thin. The absence of any sources above the 



line is consistent with both the abundance ratio values of 
Penzias (3.65) and Ladd (4.0), but not higher values. As 
so many of our sources lie close to this line it suggests 
that the value of 3.65 is the upper limit for this ratio. 
However, the majority of our sources lie below this value 
suggesting that the C 18 emission is not optically thin. 
The figure also shows the expected ratio of integrated in- 
tensities for optical depth 2 in the C 18 transition. Due 
to the overlap of the hyperfine components of the C 17 
the optical depth at the line peak depends on the velocity 
dispersion. The line on the figure has been calculated for 
a FWHM velocity of 3 kms -1 (and a C 18 to C 17 abun- 
dance ratio of 3.65) giving a corresponding C 17 peak 
optical depth of 0.5. If the FWHM velocity was 1 kms -1 , 
the C 17 optical depth would be reduced by a factor of 
0.67, whereas if the the FWHM velocity was 5 kms -1 , the 
optical depth would be increased by a factor of 1.05. The 
figure shows that the C 17 emission is not highly opti- 
cally thick towards these sources, a result consistent with 
optical depth implied by the fit to the hyperfine structure 
of the observed lines (Sec. 13. l| ). The figure also shows that 
there is no correlation between line strength and optical 
depth. 
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Fig. 2. Comparison of the integrated intensities of C ls O and 
C 17 0. The solid line denotes an abundance ratio of C lr O to 
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3.3. Column Densities 

Since the C 17 emission appears to not be optically thick, 
we calculate the column densities assuming that the C 17 
is optically thin. This is consistent w ith the analysis of 
C 17 Q towa rds similar ob jects b y both Ivan der Tak et al. 
(|2000h and lFontani et alJ (|2006h . 

To calculate the CO column densities we used the 
following general expression assuming LTE and optically 
thin emission; 



Ntat = 



3A: 



Q{T ex ) 



%TT 3 vH 2 g u exp(-E u /kT e: 



T m bdv 



Here g u is the statistical weight, Q(T) is the partition 
function and fi is the electric dipole moment for that 
molecule. Expressing the integrated intensity in Kkms -1 , 
the dipole moment in Debye (yU=0.11D for CO) and the 
frequency in GHz, this reduces to 



N tnt = 1.67x10 



14 



Q{T e . 



/i 2 (Debye)i/(GHz) 5u 



exp 



Eu 

Tea 



of < 20 K. To estimate the C 17 column density we have 
adopted 30 K, a compromise value intermediate between 
cold dust and NH3 temperatures. The column density is 
of course sensitive to the assumed excitation temperature 
and is at a minimum at T ex = 17 K. At 10 K and 30 K the 
column densities are approximately equal and about 16% 
higher than at the minimum. 

For comparison with the observations of the dust to- 
wards the sources it was necessary to re-convolve the 
850 /im data to a beamsize to match that of the CO beam. 
We first calculate the mass in the beam using the expres- 
sion 



Mi 



(beam) 



q J2 

(feeam) 



where we take the value for the opacity from WFS04 
as K85o=l-54x 10 -2 cm 2 g _1 and the dust temperature 
(Tdust) as the cold dust component given in SBSMW. We 
then derive the beam-averaged hydrogen column density 
from the submilli metre mass using the relationship from 
Hildebrandl (|1983l ) given as 



(beam) 



[ir9id 2 }N(H + H 2 )m H lJL 



T m hdv 



This method assumes a constant excitation temperature 
(T ex ) along the line of sight which is a free parameter in 
the analysis. For all their observed sources, SBSMW also 
estimated the temperature of the cold component of the 
dust; the mean over the whole sample was 45 K. SBSMW 
also measured the gas temperature towards many of the 
sources observed here using NH3 . Over their entire sample 
SBSMW find a mean temperature of 19 K and for all but 
two of the sources observed here they found temperatures 



where mn is the mass of a hydrogen atom and fi is the 
ratio of total gas mass to hydrogen mass and 0± is the 
beam radius. In both steps the distance to the source, d, 
is taken to be the near-distance if there is an ambiguity 
(SBSMW). The results for N(C 17 0) are given along with 
the hydrogen column density derived from the dust con- 
tinuum emission in Table [3j 

Figure [3] shows the column density C 17 plotted 
against that calculated for H 2 from the dust continuum 
emission. Also shown are lines of constant CO abun- 
dance chosen to constrain the data. It is immediately ev- 
ident that these graphs show a large scatter in the abun- 
dances. The C 17 data are constrained by abundances 
of ~lxlO~ 7 and ~7xl0~ 9 , a factor of approximately 
fourteen spread in the abundance. Takin g the standard 
[C 17 Q 1/[H 2 1 abundance to be 4.7xl0 -8 (jFrerking et~"ai1 
ll982T ), we would expect our sources to lie either on this 
line or below it in the case of CO depletion, yet it is clear 
that this value does not represent the upper limit for our 
sources and a significant fraction lie up to a factor of 2 
above this line. 

To explore the reliability of the derived C 17 abun- 
dances, we have also run detailed models of the emission 
expected from these sources. These models allow us to 
investigate the possible consequences of some of the as- 
sumptions made when calculating the column densities, 
particularly the assumed excitation temperature and dust 
temperature, together with the role of density gradients 
within the beam. 



4. Modelling 

To model the emission of C 17 and C 18 we used 
the ID radiative transfer code R ATRAN developed by 
Hogerheiide fc van der fail (|2000h . RATRAN utilises the 
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Fig. 3. Comparison of the column density C lr O with the 
H2+H column density derived from the dust continuum emis- 
sion. Typical uncertainties are shown in the lower right corner 
of the panel. The dashed lines represent lines of constant abun- 
dance: Short dashed = 5xl0 -7 , dot-dashed = lxlO -7 , long 
dashed = 5xl0 -8 & dotted = lxlO -8 . The circled sources 
indicate those which have been modelled (see Sec|3|. 



Monte Carlo method to calculate the radiative transfer of 
molecular lines through a dusty shell. 

IWilliams. Fuller fc Sridharanl (|2005l hereafter 
WFS05) modelled the 850 /xm emission for a selection of 
the sources observed by WFS04 using the ID radiative 
transfer code DUSTY (jlvezic. Nenkova fc Elitzurl Il999h 
in order to determine the physical parameters of the 
envelopes around the central embedded sources. Of those 
sources for which we have both C 17 and C 18 data, 14 
were successfully modelled by WFS05. We have utilised 
the parameters of the best-fit models from WFS05 to 
generate the input for RATRAN. 

WFS05 identified certain sources as asymmetric and 
used radial slices to explore the structure in various direc- 
tions, as a result four of our 14 sources actually have mul- 
tiple models in WFS05. For these sources we correspond- 
ingly generated multiple models for RATRAN. These are 
indicated by alphabetical suffixes to the source names. 

As RATRAN is a ID code it models the sources as 
spherically symmetric and hence neglects any openings 
or elongations due to outflows or winds. The model as- 
sumes a dust-free cavity surrounding the star from radius 



Table 4. Parameters used in the detailed modelling of the 
sources. Where ri is the inner radius of the modelled envelope, 
T2 is the outer radius and a is the power-law index of the 
number density profile n(r) where n oc r~ a . Column 5 gives 
the hydrogen number density at r 1 . 



r = to r = T\. Between ri and r 2 lies the envelope 
to be modelled. We divided the radial distance between 
ri and r 2 into 30 logarithmically spaced shells to give in- 
creased resolution towards the inner region of the shell 
where the temperature and density profiles change most 
rapidly. Thi rty shells were used corres ponding to the num- 
ber used bv lvan der Tak et all (|200fj| . hereafter WEB) in 
their RATRAN modelling of massive young stars. We also 
ran a number of identical models using both 15 and 60 
shells and found a negligible effect on the output com- 
pared to the 30 shell models. 

To determine the input for RATRAN we used the fol- 
lowing parameters derived in WFS05: power-law index (a) 
of the envelope density profile n(r) where n oc r~ a , the 
temperature at the inner boundary T%, the scale of the 
envelope defined by the inner radius, r\ and Y, where Y 
is the ratio between r 2 and r\ . The outer radius can then 
be calculated from r 2 = Yt\. 

We generated the density profile for each model by 
combining the power-law index (a) of the envelope den- 
sity profile given in WFS05 with the density at r% which 
we calculate from the 850 /im mass and volume of the en- 
velope (WFS04). As the dust temperature profile is not 
given in WFS05 we regenerated the WFS05 models using 
DUSTY to determine the individual temperature profile 
for each model. This profile was then interpolated to the 
radii of the 30 shells used to trace the envelope. Given the 
high densities in these regions we have assumed that the 
gas kinetic temperatures and the dust temperatures are 
equal. 

The turbulent line widths were taken from Gaussian 
fits of the optically thin C 17 data. The presence of hyper- 
fine structure for the C 17 was neglected for the models. 
By fixing the linewidth from our observations we were able 
to isolate the CO abundance relative to hydrogen as the 
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Fig. 4. Comparison of the C 17 J = 2 — + 1 (upper spectra) and C ls O J = 2 — » 1 (lower spectra) spectra with the results of 
the best fit results of the RATRAN models (solid curves). The model parameters are given in Tables [4] and [S] 



free parameter. For simplicity we restricted our models to 
a simple uniform abundance of C 17 throughout the enve- 
lope of each source. The input parameters for each model 
are given in Tabled] 

These parameters were then used in co njunction with 
the ado pted collisional rate coefficients of iFlowerl (|200lh 
(see also Schoier et al . 2005) to generate the predicted line 
emission. The C 17 abundances were adjusted to produce 
the best match to the observed C 17 line profiles. The 
modelled sources are circled in Figure [3l illustrating that 
they cover a range dust and gas column densities and a 
range of implied C 17 abundances. 

The output spectra from RATRAN were continuum 
subtracted and convolved with a Gaussian with a FWHM 
equal to the appropriate JCMT beam for comparison with 
the observations. We analysed the quality of the models by 



using the reduced x 2 method, aligning the central channel 
of the model spectra with the peak of the Gaussian fit for 
the data. Figure [4] compares the observations with the 
results of the best fit models. 



4.1. Modelling Results 

Figure O shows the success with which these (simple) mod- 
els can match the observed C 17 line profiles and also the 
sensitivity of the fit to the assumed C 17 abundance. In 
this case the C 17 abundance is constrained within the 
noise level of the observations to within ~ 10%. 

Since RATRAN calculates the continuum emission in 
parallel with the line emission it was possible to check on 
the consistency of the models by comparing the predicted 
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Fig. 5. The observed C 17 J = 2 -» 1 emission towards 
WFS14 overlaid with the results of three models. The only 
difference between the three models is the assumed C 17 abun- 
dances of 5xl0 _s , 5.5xl0 -8 & 6xl0 -8 . The abundance is con- 
strained within the noise level of the observations to within 
-10%. 



850 /im flux with that measured. This is especially impor- 
tant in this work where we are looking at the ratio of the 
line to continuum emission. For a number of sources mod- 
elled the predicted and measured 850 flux differed by 
more than a factor of 3. This is not completely surprising 
as the best fit models of WFS05 were selected to provide 
the best overall fit to the source SED and 850 /im distri- 
bution, as such they do necessarily provide excellent fits 
to the 850 /im flux alone. Therefore we have only consid- 
ered the 8 sources where the predicted 850 /im continuum 
emission differs from the measured value by less than a 
factor of 3. The results of these sources are given in Table 

El 

4.1.1. C 17 Analysis 

A comparison between the observed line profiles and the 
best fit model results for all the modelled sources is 
shown in Figure [4] The models can successfully match 
the observed C 17 line profile with a reasonable range of 
C 17 fractional abundances ranging over the sample from 
1.8xl0~ 8 to 1.5xl0~ 7 . Importantly for the analysis of the 
sources for which models of the continuum emission do not 
exist, comparison of the C 17 abundance inferred by the 
modelling and that derived directly from the observations 
agree within a factor of —2 (Table [5]). This result indi- 
cates that the assumptions made in directly estimating 
the C O abundance from the observations are not overly 
biasing the results. 

4.1.2. C 18 Analysis 

For completeness we also attempted to model the C 18 
lines for the sources modelled in C 17 0. Initial modelling 



of the C 18 lines involved taking the best-fit C O abun- 
dance and scaling it up by a range of potential values of 
the C 18 to C 17 abundance ratio; these were R=2.8, 
3.65, 4.0 & 4.5 (as mentioned in Sec. 13. 2\i . All other free 
parameters remain constant with the exception of the ve- 
locity dispersion which was adopted from the C 17 hfs 
fit (as the best estimate of the intrinsic velocity disper- 
sion within these regions). The results are also shown in 
Figure 2J On inspection it is clear that a number of the 
models suffer from high optical depths, resulting in a flat- 
tening and broadening of the line shape, which is not seen 
in the observations. The inconsistencies between the data 
and the models for this isotopologue is perhaps not unex- 
pected given the two velocity component profiles towards 
many of the sources which suggest that more than a simple 
envelope is contributing to the line profile. 

An example of such a scenario would be a central core 
(as represented in our models), but surrounded by an ex- 
tended low density envelope with the C 18 having differ- 
ent optical depths in the envelope and in the core. The 
success of modelling the C 17 from the core (Sec. 14. ip 
would suggest that the core has the higher optical depth. 

For sources with two component C 18 line profiles 
(Table [7]), the C 17 line is often intermediate in width 
between the two C 18 components. This could also point 
towards the presence of two different components of ma- 
terial. If there are indeed two components then the simple 
column density analysis and model are actually overesti- 
mating the actual C 17 and C 18 column densities in the 
core, as some of the column density attributed to the core 
in a single component model is actually associated with 
the second component. This suggests that the abundances 
of these species in the core could be even lower than the 
values derived here indicate. 

It is possible to test this prediction with maps of these 
sources; C ls O emission extending beyond the detected 
850 /im emission would provide evidence for an outer enve- 
lope. It is possible that the relatively small chopper throw 
of the 850 /im observations may have artificially removed 
this component from the SCUBA observations. For some 
sources this resulting underestimation of dust emission 
could account for the high abundances seen towards a few 
sources. 



5. Discussion 

The determination of the C 17 column densities and 
hence abundances contain a number of sources of uncer- 
tainty. Since both the C 17 lines are detected with good 
signal to noise ratios (TableHJ Figured]), and the dust con- 
tinuum fluxes of the sources are well measured (WFS04) , 
the statistical uncertainties due to noise are small. The 
calibration uncertainties may be larger, but are likely to 
affect all the sources equally. However the results do also 
depend on a number of other assumptions. Ideally all the 
sources observed should be modelled in detail to deter- 
mine their C O abundance throughout their envelopes 
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Table 5. Parameters of the best fit RATRAN models. 
Column 2 is the fractional abundance of C lr O required for the 
model to match the observed J = 2 — > 1 line profile. Column 3 
is the reduced x 2 between the observed C 17 J = 2 — > 1 tran- 
sition and the modelled emission. R(H2) is the ratio of mea- 
sured 850 /im flux to that generated by RATRAN. R(C ir O) 
is the ratio of the best fit abundance from the modelling to 
the abundance derived from comparing the data to the 850 /J,m 
hydrogen abundance: [C 17 0/H2] mo d/[C 17 0/H2]data. Column 4 
presents the best fit ratio of C ls O/C 17 needed to match the 
C 18 observations, scaling up the best-fit C 17 abundances 
listed. The model line profiles are compared with the observa- 
tions in Figure [4] 



and cores. However temperature and density profiles are 
not (yet) available for all the objects, so this is not cur- 
rently practical. Nevertheless the objects which have been 
modelled in detail do suggest that systematic issues do 
not dominate the results. The similar C 17 abundances 
derived from the modelling of a subsample of the sources 
which span a range of source properties and the C 17 
abundances derived directly from the observations, sug- 
gests that despite the approximations and assumptions 
it requires, the direct derivation produces reasonable esti- 
mates of the average C 17 abundance towards the sources. 
It is these abundances on which we will focus. 

To determine the origin of the scatter seen in Figure 
[3] we must investigate the possible causes ranging from 
underlying physical differences to the simplifying assump- 
tions. The consistency within a factor of 2 of the C 17 
abundances generated using RATRAN and those derived 
straight from the data suggests that the scatter between 
the observed C 17 column density and H 2 column density 
is not a result of an overly simplistic data analysis or the 
assumed excitation temperature. 

It is important to recognise that the hydrogen column 
density derived from the dust emission suffers from uncer- 
tainties due to the assumed dust properties. In calculating 
the dust column density we have adopted a constant value 
of k, the dust mass opacity, for all t he sources. For our cal- 
culati ons we adopted a value from lOssenkopf &: Hennine 
( 19941 ) assuming a gas density of 10 6 cm -3 , with thin 
ice mantles and 10 5 years of coagulation. However we 
investigated the influence of this assumption by recalcu- 
lating N(H2) for the full range of possible opacities from 



Ossenkopf & Hcnning. By manipulating the opacity it is 
possible to move the entire sample in Figure [3] to a range 
either above the canonical C 17 abundance of 4.7xl0 -8 
or below it. However assuming that all the sources should 
lie at this value, it is clear that changing the dust in or- 
der to allow that to happen would have to be done on a 
source by source basis and utilising the full range of po- 
tential opacities. 

For those sources that lie to the left of the canoni- 
cal abundance line, an opacity corresponding to grains 
that have not undergone any coagulation would need to 
be adopted. This lack of grain coagulation would seem to 
imply that these sources, and the dust surrounding them, 
are young, not yet having had time for the dust to coag- 
ulate. To the other side of the line lie the sources which 
appear to be the most depleted of our sample. To change 
their abundances would require an opacity corresponding 
to some degree of coagulation and for the most under- 
abundant sources this would need to include an absence 
of ice mantles. The scenario of coagulated grains without 
ice would be consistent with the dust close to a star which 
has heated its surroundings and evaporated all icy man- 
tles. However the temperature profiles derived by DUSTY 
modelling show that only a small part of the envelopes arc 
this warm, with the bulk of each envelope lying at tem- 
peratures below 50 K. Under these conditions we would 
consider the complete absence of ice mantles to be an un- 
reasonable assumption. The dust models indicate if man- 
tles of any kind are present, they would limit the change 
in the inferred hydrogen column density to not more than 
12% of the calculated values. 

In addition, van der Tak et ahl ( 19991 ) found that by 
comparing the dust mass from FIR and submm data to 
the gas mass derived from regions without depletion of 
C 17 0, only the opacities for dust grains with ice man- 
tles reproduced the standard dust-to-gas ratio of 1:100. 
Additionally figure shows a number of the sources hav- 
ing abundances in excess of 4.7x 10 -8 , which if not due to 
particular circumstances for these objects, could indicate 
that for the sample as a whole the hydrogen column den- 
sity could be underestimated. Given these points we con- 
clude that the scatter seen between the abundances does 
indeed reflect physical differences between the sources. 

In principle selective photodestruction of the C 17 in 
the inner regions of the sources close to the central heating 
sources could give rise to the varying C 17 abundance. 
However selective photodestruction can only be an issue 
in the region up to optical depth ~ 1 in the UV from the 
inner edge of the circumstellar envelope. Although this 
optical depth can be contributed by either dust or line 
self-shielding, since the circumstellar envelopes are dusty, 
it is the dust that will be the dominant UV opacity source. 
This UV opacity corresponds to a visual optical depth 
~ 0.1 in the optical. Since this corresponds to about 1% 
of the thickness of the envelope (WFS04), the internal UV 
flux can not be affecting the C 17 abundance by more 
than - 1%. 
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Table 6. Mass-weighted temperatures in Kelvin, calculated 
for those sources for which the data are available. 



The low optical depth of the C 17 transitions sug- 
gest that if C 17 is uniformly abundant throughout each 
source, the C 17 emission, and hence C 17 column den- 
sity should, like the dust emission, trace the total column 
density even though the C 17 has a relatively low critical 
density. 

The most likely explanation for the variation in the 
C 17 abundance appears to be that it reflects different 
gas phase abundance of C 17 due to its freeze out onto 
dust grains. Such depletion of C 17 (and by implica- 
tion, CO and possibly other species) towards similar high 
mass sources has be en previously reported by WEB and 
Fontani et all (|2006l ). 



5.1. Depletion of C" O 

If the scatter in the C 17 abundance is indeed due to 
depletion, then one might expect that the regions with 
the coolest dust would have the highest depletion. WEB 
investigated this using the mass-weighted temperature, 
T m w, to define a characteristic temperature for each 
source. 

Figure [5] shows the C 17 abundance against the mass- 
weighted temperature for all the sources for which the 
necessary parameters were available. The mass- weighted 
temperature is defined by WEB and was calculated here 
using the temperature profile from the DUSTY modelling 
of the sources. The figure shows data points from both 
our sources and those from WEB. The two data sets are 
in reasonable agreement, with the distribution of both be- 
ing consistent with depletion at low temperatures. The 
Spearman correlation coefficient indicates that the C 17 
abundance and mass- weighted temperature are correlated 
at a 96% confidence level. 

We do however see a number of our sources exhibiting 
higher abundances at lower temperatures than WEB, but 
none with showing the reverse. This is consistent with an 
absence of depletion at higher temperatures, whilst the 
sources at low temperature and high abundance could be 
explained by having additional contributions of C 17 for 
regions outside the modelled core. This would be consis- 
tent with the C 17 linewidths lying intermediate between 
the two components of the C 18 lines as discussed above. 

The degree of depiction in a region can provide an es- 
timate of the time for which gas has been cold and at 
high density. This provides a measure of timescales for 




Fig. 6. C ir O (J = 2 — » 1) abundance against the mass- 
weighted temperature, T mw . The solid dots denote our data 
whilst the circles represent the data from Figure 11 in WEB. 



the formation of high-mass stars. Models of the phys- 
ical and chemical changes occurring in the gas and on 
dust grains during the pre-protostellar phase of the evo- 
lution of cores, lead to predictions for the expected abun- 
dances of molecules such as CO as a fu nction of time 
dViti fc WiUiamslll999l iFlower et al. I l2005l ). With detailed 
knowledge of depletion levels they become a viable tool 
for determining the ages of the cores. 

The detection of both C 17 and C 18 towards all of 
the sample shows that that there are no sources towards 
which essentially all CO has been depleted. This places 
an upper limit on the time for which the dust can have 
existed at low temperatures of T < 10 K. As we do not 
find depletion to exceed more than a factor of ~10 we can 
use the abundances predicted by IFlower et al] (|2005h to 
estimate that the cores being traced are not older than 
^3.5xl0 5 years. (For a depletion factor of 5 this reduces 
to ~2.2xl0 5 years). 



5.2. Evidence of Source Evolution? 

This set of source s has also been observed by 



Beuther et al.l (|2002af ) at 1.2 mm continuum emission and 
a number of CS lines. Figure [7] shows a comparison of 
the linewidth of C 17 J = 2 — > 1 (determined from the 
fit to the hyperfine structure of the transition) against 
the linewidths of CS J = 3 -> 2 and C 34 S J = 3 -> 2. 
It is clear that at smaller linewidths the two are most 
closely correlated, with approximately equal linewidths. 
However the linewidths tend to have a larger ratio at larger 
linewidths. For CS linewidths < 3 kms -1 the ratio of the 
CS linewidth to the C 17 linewidth is 1.17 ± 0.05, while 
for CS linewidths > 3 kms -1 , the ratio is 1.56 ±0.09. The 
same trend is also seen in the rarer species C 34 S, where 
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Fig. 7. Plot showing the C lr O hfs linewidth against (a) 
the C 34 S (J = 3 -» 2) and (b) the CS (J = 3 -> 2) linewidths 
jBeuther et alj|2002af ). The solid lines show the unity ratio. 



the ratio for C 34 S linewidths < 3 kms" 1 is 1.00 ± 0.05 
and 1.22 ± 0.07 for the larger linewidth sources, suggest- 
ing this is not a result of high optical depth in the CS 
line. We speculate that this increasing ratio of CS to C 17 
linewidth could arise as a result of CS being more affected 
than the C 17 by heating and stirring of the material by 
the central source, through the action of the outflow. 

Interestingly there appears to be a connection be- 
tween the C 17 abundance towards a source and its CS 
linewidth. Comparing the distributions of C 17 abun- 
dance for those sources with Aw(CS) > 3 km s _1 and those 
for which At>(CS) < 3kms -1 , a Kolmogorov-Smirnov two- 
sample test indicates a 98.6% likelihood of a difference 
between the two groups (see Figure UK). We also con- 
sidered just the upper (Aw(CS) > Skms" 1 ) - 30% of 
the sample, comparing its distribution of C 17 abundance 
with that of the lower (Au(CS) < Skins" 1 ) - 30% of 
the sample. Repeating the K-S test for these groups of 
sources increases the statistical significant of the differ- 
ence between these new groups to 99.7% (see Figure [8b), 
despite the smaller number of sources in each group. This 
increase in the difference between the groups might be 
expected if these groups are not distinct categories but 
rather represent two extremes in a continuum of source 
properties, from small linewidth, high C 17 abundance 
objects to large linewidth, low C 17 abundance objects, 
perhaps tracing an evolutionary progression. 

To further probe the statistical significance of this dif- 
ference between the sources we performed two tests. In the 
first we randomly shuffled the C 17 abundances amongst 
the sources. The sample was then again divided between 
objects with Aw(CS) > 3kms _1 and those for which 
Aw(CS) < 3kms _1 . A K-S test was then performed com- 
paring the resulting distributions of C 17 abundance for 
the two linewidth groups. Repeating this randomisation 
and testing 100 times we found no cases where the sta- 
tistical difference between the shuffled abundances exceed 
that of actual data. In other words the difference in C 17 
abundance between the sources with large linewidths and 
those with small linewidths has a probability of less than 
1% of being a chance coincidence. 



For the second test we take into account the obser- 
vational uncertainties associated with the abundances by 
randomly resampling the abundances associated with each 
source. For each source we regenerated a C 17 abundance 
by randomly drawing an abundance from a Gaussian dis- 
tribution with a dispersion equal to the uncertainty in the 
measured C 17 abundance of that source. We did this for 
the whole sample, split the sample on the basis of the CS 
linewidth and again apply a K-S test to intercompare the 
resulting distributions of C 17 abundance. This was re- 
peated 100 times. We found that on average the two sam- 
ples differed at 98% confidence level, with the confidence 
level never less than 97%. This indicates that the differ- 
ence in abundance between our two groups is unlikely to 
arise as a result of the statistical uncertainties. 

The association of low C 17 abundance with sources 
with large linewidths is somewhat surprising. However it 
can be understood if the low C 17 abundance (which is 
calculated integrated across the whole line profile and av- 
eraged across the telescope beam) reflects the conditions 
in the cold environment in the outer, extended envelope 
around the central source, while the broad linewidth re- 
flects some component close to the central source, possibly 
related to outflow. Mapping the spatial distribution of the 
depletion towards these objects by comparing molecular 
line and deep dust continuum maps would be able to test 
this interpretation. 

Whilst the central regions of these regions certainly 
have temperatures reflecting the presence of a luminous 
embedded source, the models of WFS05 imply that the 
majority of the extended envelope lies at temperatures 
below 50 K. In many cases the models show the tempera- 
ture dropping below 20 K towards the outer regions of the 
envelopes, temperatures which are similar to those derived 
from the NH3 observations of SBSMW. It is presumably 
in these outer regions, far from the heating and effects 
of the outflow from the central source, that the C 17 is 
frozen out onto the surface of dust grains. 

The high C 17 abundances displayed by those sources 
with the lowest linewidths could be consistent with two 
possible scenarios; either these sources are extremely 
young and depletion has not yet occurred to a significant 
degree, or alternatively the sources are much older having 
undergone substantial heating by an embedded population 
of sources, or shocks from outflows, which has evaporated 
the ice mantles containing the C 17 back into the gas 
phase. 

This second scenario could be confirmed by the posi- 
tive detection towards these sources of the molecules which 
trace hot core emission. SBSMW have searched many of 
these sources for such species and shown that at least some 
of these sources are evolved enough to have locally heated 
their natal cores to temperatures of >100K. However the 
failure to detect these molecules towards some sources im- 
plies that these sources are either highly evolved with the 
hot core molecules have been destroyed by ion reactions 
and ages ~10 5 years (e.g. Hatchell et al. 1998), or else 
would suggest that these are extremely young sources. A 
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Fig. 8. Distribution of C 17 abundance, (a) The solid bars 
show the distribution of abundances for those sources with 
Aw(CS) > 3 kms" 1 whilst the hollow bars represent those 
sources with Av(CS) < 3kms _1 . (b) The solid bars show 
the distribution of abundances for the top 30% of the sam- 
ple (Aw(CS) ^> 3 kms -1 ) while the hollow bars show the lower 
30% of the sample (Aw(CS) < 3 kms -1 ). 



more det ailed discussion of these possible scenarios is pre- 
sented in lThomas fc Fullerl (|2007l ). High spatial resolution 
observations of these regions would help to distinguish be- 
tween the two pictures by identifying the location of the 
depleted material. Very young cores would be expected 
to be depleted in C 17 in their centres where the den- 
sity is highest but to have normal C 17 abundance in 
their more extended, lower density envelopes. For a more 
evolved core which has formed a high mass protostar at 
its centre, this central object will be heating the interior 
region of the core, releasing frozen out C 17 back into the 
gas phase. However sufficient time may have passed for 
C 17 to be depleted in the outer parts of the core. 



6. Summary 

We have presented spectra of 84 candidate HMPOs in the 
J = 2 — ► 1 line of C 17 and for a selection of these we have 
also observed C 18 J = 2 — > 1 and C 17 J = 3 -> 2. We 
have used the data to calculate the C 17 column densities 
and for those sources for which models of dust emission ex- 
ist we constructed models of the CO emission to check the 
validity of our assumptions regarding excitation tempera- 
ture and dust temperature. We have found the following: 

The ratio of C 17 column density to H2 derived from 
850 /im dust emission has a significant scatter across our 
sample. The scatter between the most extreme sources is 
~14 and is consistent with the abundance of C 17 differ- 
ing on a source- by-source basis. 

Our derivation of the H2 column density assumed a 
constant value for k for all the sources, however there are 
no corrections to k which when applied globally to the 
sample to eliminate this scatter. Indeed the elimination of 
the scatter would require the sources lying at the extremes 
of the scatter to possess dust properties which are physical 
unlikely. 



The method and assumptions used to calculate the 
C 17 column densities were tested by modelling a sub- 
set of these sources with the radiative transfer code 
RATRAN. For all suitable models the best-fit abundance 
of C 17 matched the data to within a factor of ~2, sup- 
porting our assumptions and removing the possibility of 
the assumptions accounting for the range of abundances. 
We conclude that this scatter arises as a result of deple- 
tion present in a number of these sources. However the 
maximum degree of depletion we find does not exceed 
~7. This range of depletion suggests that outer regions 
of these sources have lifetimes during which they are cold 
and dense of about 2 — 4 x 10 5 years. 

We derive the mass-weighted temperature and find a 
correlation with C 17 abundance, consistent with WEB 
who from this correlation also infer the presence of deple- 
tion towards sources similar to those discu ssed here. 

The se sources have been studied by iBeuther et al. 
(|2002al ) and a comparison of our C 17 linewidths to 
their CS J = 3 — > 2 linewidths shows a division between 
those which have Aw(CS) > 3 kms -1 and those for which 
Aw(CS) < 3 kms -1 . A statistical comparison of the C 17 
abundances of these two groups reveal a trend of the 
first group displaying significantly lower abundances than 
those in the second group with a statistical significance 
exceeding ^98%. 

We suggest that this range of C 17 abundance and the 
implied degree of depletion between the sources may re- 
flect a spread in evolutionary status amongst the sources. 
Higher angular resolution observations of dust emission, 
C 17 and other species towards these objects will be 
important in confirming this suggestion. It is possible 
that high resolution observations will identify regions with 
higher degrees of depletion than can be measured with the 
relatively large telescope beam used for the observations 
presented here. 

Depletion towards these objects shows that during the 
evolution of these cores the gas has remained cold and 
dense for long enough for the trace species to deplete. More 
detailed models of how trace molecules deplete in massive 
dense cores, combined with higher angular resolution ob- 
servations, should be able to provide tighter constraints 
on the lifetime of these regions, including how long they 
exist before they form massive stars. 
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Table 2. Observed line parameters. The columns give the peak intensity, velocity, linewidth and integrated intensity for C 17 0. 
The first three parameters are derived from hyperfine fitting of the data using CLASS. The integrated intensity was obtained 
using a single Gaussian fit. Three sigma upper limits are shown for sources where a line was not detected above this threshold. 
An * indicates an data that has been been observed in the previous year and is given in Table [T] 
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Table 3. Column densities calculated from the line intensities with an average rms of ~0.2K. The hydrogen column density 
corresponds to a 20" beam to directly compare with the J = 2 — > 1 transitions of each isotope. 
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2.43 


60.40 


2.46 


58.63 


17 


3.02 


45.12 


2.65 


45.09 


2.72 


45.08 


2.20 


44.87 


3.70 


45.67 


20 


2.99 


34.35 


2.70 


34.29 


2.82 


34.39 


1.60 


34.99 


2.97 


34.03 


21 


2.76 


34.40 


2.25 


34.35 


2.25 


34.39 


1.97 


34.46 


3.94 


33.61 


29 


2.59 


59.28 


2.18 


59.28 


2.69 


58.94 


1.83 


58.42 


1.93 


59.93 


30 


3.05 


95.64 


2.96 


95.75 


2.44 


95.69 


1.46 


95.65 


5.80 


95.79 


39 


3.08 


96.07 


2.61 


96.03 


2.64 


96.07 


1.55 


96.36 


3.72 


95.53 


78 


2.82 


21.84 


2.47 


21.80 


2.83 


21.72 


0.84 


21.70 


3.10 


21.73 


79 


2.91 


22.60 


2.33 


22.54 


2.85 


22.53 


1.30 


22.49 


3.96 


22.66 


90 


3.25 


-3.63 


2.94 


-3.67 


2.83 


-3.64 


1.28 


-3.89 


3.59 


-3.41 


108 


3.06 


-53.06 


2.87 


-53.10 


2.79 


-53.18 


1.50 


-53.06 


4.16 


-53.37 


109 


2.56 


-44.35 


2.33 


-44.38 


2.74 


-44.29 


1.80 


-44.58 


3.57 


-43.82 


110 


2.97 


-55.00 


2.09 


-54.74 


2.31 


-54.71 


1.35 


-54.66 


3.41 


-54.24 


112 


2.03 


-18.53 


1.55 


-18.62 


1.84 


-18.45 


1.18 


-17.41 


1.31 


-18.65 



Table 7. All columns (except column 1) in units of kms -1 . The Gaussian linewidths were obtained using SPECX, whilst the 
hyperfine linewidths were obtained using the hfs method in CLASS. Multiple Gaussian components are only given for those 
sources where appropriate. 



H. S. Thomas and G. A. Fuller: HMPOs: Depletion 



WFS 


IRAS NAME 




R.A. 


Dec. 




DATES OF OBSERVATIONS 










(J2000) 


(J2000) 


C ls O(2-l) 


C 17 0(2-l) 


C 17 0(3-2) 


WFS1 


IRAS 


05358+3543 


05 


39 


10.8 


+35 45 


16 




05/07/91 

UO / U 1 1 £i -L 


05/07/23 


WFS3 


IRAS 


05490+2658 


05 


52 


11.0 


+27 00 


34 




05/07/21 

W ' J 1 W 1 / — 1 _L 


05/07/23 


WFS4 


IRAS 


05490+2658 

V / O 1 • J \J | — \ J O W 


05 


52 


12.1 


+27 00 


11 




05/07/21 

w *-J / w t / £j J. 


05/07/23 


WFS6 


IRAS 


05553+1631 


05 


58 


13.4 


+16 32 


00 




05/07/21 

W * J 1 w i j i 


05/07/23 


WFS7 


IRAS 


18089-1732 


18 


11 


45.2 


-17 30 


43 






05/07/07 


WFS8 


IRAS 


18089-1732 


18 


11 


51.5 


-17 31 


34 






05/07/07 


WFS 11 


IRAS 


18089-1732 


18 


11 


57.0 


-17 29 


34 






05/07/07 


WFS12 


IRAS 


18090-1832 


18 


12 


02.1 


-18 31 


58 




05/08/03 




WFS13 


IRAS 


1 81 02-1 800 

lO lUi lOWW 


18 


13 


11.7 


-18 00 


04 




05/08/03 

UO 1 UO / UO 




WFS14 


IRAS 


181 51-1208 

lO lO 1 IZiWO 


18 


17 


58.2 


-12 07 


28 


04/05/1 3 

U4: / UO j -LO 


04/05/13 04/05/14 

U4:/UU/ lO; U4:/UU/ 14: 




WFS15 


IRAS 


18159-1550 


18 


18 


48.4 


-15 49 


00 




05 /08 /03 

WO 1 WW/ UO 




WFS16 


IRAS 


1 81 89_1 4S3 

lO lOi 14:00 


18 


21 


08.9 


-14 31 


46 


04/05/1 3 

U4: / UU / lO 


04/05/1 3 

U4: / UO j -L O 




WFS17 


IRAS 


18223-1243 


18 


25 


10.6 


-12 42 


27 


04/05/13 

v. / i j wo/ -LO 


04/05/14 

04: / \J ' J j i i 




WFS18 


IRAS 


18247-1147 


18 


27 


31.4 


-11 45 


55 




05/07/27 05/08/03 




WFS19 


IRAS 


18264-1152 


18 


29 


14.3 


-11 50 


22 




05/07/27 05/08/03 

WO / W J / J - wo / uu / wo 




WFS20 


IRAS 


18272-1217 

_L W w I Zj 1L1 i 


18 


30 


02.2 


-12 15 


40 


04/05/1 3 

VJrt f WO/ -LO 


04/05/14 

Wt: / WW / It: 




WFS21 


IRAS 


1 8272-1 21 7 


18 


30 


03.2 


-12 15 


11 


04/05/1 3 

U4: / UU / lO 


04/05/14 

U4 / WW / -L4: 




WFS22 


IRAS 


18290-0924 


18 


31 


43.4 


-09 22 


26 


04/05/14 


04/05/14 

W4: / WW/ i i 


05/07/22 


WFS23 


IRAS 


1 8290-0924 


18 


31 


44.0 


-09 22 


17 




05/08/03 

WW / WW / wo 




WFS 24 


IRAS 


1 8306-0835 

iwoww wwoo 


18 


33 


17.3 


-08 33 


28 




05/08/03 

WW / WW / wo 




WFS25 


IRAS 


18306-0835 


18 


33 


23.9 


-08 33 


33 


04/05/13 

V. / I / WO/ _LO 


04/05/13 

W i / WW/ -L w 




WFS28 


TR AS 


luOlU Uu^J 


18 


33 


47.9 


-08 23 


52 




UO / UO / WO 




WFS29 


TR AS 


lOOO i U i 4:0 


18 


36 


27.9 


-07 40 


25 


04/0^/1 4 


04/0^/14 

U1 / UO / 11 


05/07/22 


WFS30 


IRAS 


18345-0641 


18 


37 


16.8 


-06 38 


35 


04/05/13 


04/05/13 




WFS33 


TR AS 


1 8^48 flfil (S 

-L OOIO UU1U 


18 


37 


30.5 


-06 14 


13 




0^/07/97 n^/ns/n^ 




WFS 34 


IRAS 


1 8S72-0541 

-LOO 1 WutI 


18 


37 


16.8 


-06 38 


35 


04/05/14 04/05/15 

WT: / WW/ It:^ Wt: / WO I 1 <J 


04/05/14 

Wt: / WW / It: 




WFS35 


IRAS 


18385-0512 


18 


41 


12.8 


-05 08 


58 




05 /08 /03 

WW/ WW/ WW 




WFS36 


IRAS 


1 8426-0204 


18 


45 


12.1 


-02 01 


10 


04/05/14 04/05/15 

WT: / WW/ It:^ Wt: / WO I 1W 


04/05/14 

Wt: / WW / It: 




WFS37 


IRAS 


18431-0312 


18 


45 


45.5 


-03 09 


21 




05/07/97 05/08/03 

UO / Ul 1 Zi I ^ uo/uo/uo 




WFS38 


IRAS 


18437-0216 

1Ut:'J 1 WZ-1W 


18 


46 


21.8 


-02 12 


20 




05/07/27 05/08/03 

WO / Wl / ZjI . WO / WW / wo 




WFS39 


IRAS 


18437-0216 


18 


46 


22.4 


-02 14 


16 


04/05/14 04/05/15 

WT: / WW/ It:^ Wt: / WO j 1 O 


04/05/14 04/05/15 

W4:/ WO/ 14) 04/ WO/ -LO 


05/07/07 


WFS41 


IRAS 


18440-0148 

1W4:T:W W 1 1 V. ' 


18 


46 


33.3 


-01 44 


52 




05/08/03 

WW / WW / wo 




WFS42 


IRAS 


18440-0148 


18 


46 


36.5 


-01 45 


22 




05 /08 /03 

WW/ WW/ wo 




WFS51 


IRAS 


18460-0307 

1 UtiUU WOW i 


18 


48 


37.8 


-03 03 


48 


04/05/14 04/05/15 

WT: / WW/ It:^ Wt: j WO j 1 O 


04/05/14 04/05/15 

W4:/ WO/ 14) 04/ WO/ lO 


05/07/07 


WFS55 


IRAS 


1 8472-0022 


18 


49 


52.4 


-00 18 


59 




05/08/03 

UO / UO / UO 




WFS57 


IRAS 


18488+0000 


18 


51 


24.4 


+00 04 


39 




05/07/26 

WW / W 1 / 




WFS58 


IRAS 


18488+0000 

10400 | W W W 


18 


51 


25.5 


+00 04 


11 




05/07/26 

WW / W I / iJ w 




WFS59 


IRAS 


1 8521 +01 34 


18 


54 


40.6 


+01 38 


05 




05/08/03 

Uw / UO / UO 




WFS60 


IRAS 


18521+0134 

1 ( ) o _ 1 | V / 1 *7 1 


18 


54 


44.4 


+01 37 


00 




05/08/03 

WW / WW / wo 




WFS61 


IRAS 


1 8530+021 5 


18 


55 


33.7 


+02 19 


09 




05/08/03 

UO / UO / UO 




WFS62 


IRAS 


1 8540+0220 


18 


56 


36.6 


+02 24 


45 




05/07/30 

WW 1 \J t 1 WW 




WFS63 


IRAS 


18540+0220 


18 


56 


40.1 


+02 25 


30 




05/07/30 

WW / W 1 / WW 




WFS64 


IRAS 


18553+0414 


18 


57 


53.5 


+04 18 


16 




05/07/30 

WW / W 1 / WW 




WFS66 


IRAS 


1 qni 2+0536 


19 


03 


45.3 


+05 40 


43 




05/08/03 

UO / UO / UO 




WFS67 


IRAS 


19035+0641 


19 


06 


01.5 


+06 46 


35 




05/08/03 




WFS68 


IRAS 


19074+0752 


19 


09 


53.4 


+07 57 


12 




05/07/29 




WFS69 


IRAS 


19074+0752 


19 


09 


53.9 


+07 56 


55 




05/07/29 




WFS 70 


IRAS 


19175+1357 


19 


19 


48.6 


+14 02 


26 




05/07/27 




WFS71 


IRAS 


19175+1357 


19 


19 


48.8 


+14 02 


46 




05/07/26 




WFS 72 


IRAS 


19217+1651 


19 


23 


58.6 


+16 57 


38 




05/07/27 




WFS 74 


IRAS 


19266+1745 


19 


28 


55.5 


+17 52 


00 




05/07/30 




WFS 75 


IRAS 


19282+1814 


19 


30 


23.1 


+18 20 


22 




05/07/30 




WFS 76 


IRAS 


19282+1814 


19 


30 


29.7 


+18 20 


37 




05/07/30 




WFS 77 


IRAS 


19403+2258 


19 


42 


28.8 


+23 05 


03 




05/07/28 


05/07/07 


WFS 78 


IRAS 


19410+2336 


19 


43 


10.6 


+23 45 


02 


04/05/14 


04/05/14 


05/07/07 



continued on next page 
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WFS IRAS NAME 


R.A. 


Dec. 


DATES OF OBSERVATIONS 




(J2000) 


(J2000) 


C 18 0(2-l) C 17 0(2-l) C 17 0(3-2) 



WFS79 


IRAS 


19410+2336 


19 


43 


11 


2 


+23 


44 


06 


04/05/14 


04/05/14, 04/05/23 


05/07/07 


WFS80 


IRAS 


19411+2306 


19 


43 


17 


6 


+23 


13 


57 




05/07/28 


05/07/07 


WFS81 


IRAS 


19413+2332 


19 


43 


26 


3 


+23 


40 


26 




05/07/28 


05/07/07 


WFS82 


IRAS 


19413+2332 


19 


43 


29 





+23 


40 


19 




05/07/29 


05/07/07 


WFS83 


IRAS 


19471+2641 


19 


49 


10 


1 


+26 


49 


10 




05/07/29 


05/07/07 


WFS84 


IRAS 


19471+2641 


19 


49 


11 


8 


+26 


49 


38 




05/07/29 


05/07/07 


WFS85 


IRAS 


20051+3435 


20 


07 


04 


5 


+34 


44 


45 




04/05/23 




WFS86 


IRAS 


20081+2720 


20 


10 


12 


6 


+27 


29 


13 




05/07/26 




WFS87 


IRAS 


20081+2720 


20 


10 


13 


3 


+27 


28 


21 




04/05/23 




WFS88 


IRAS 


20081+2720 


20 


10 


16 





+27 


28 


12 




04/05/23 




WFS89 


IRAS 


20081+2720 


20 


10 


18 


7 


+27 


27 


18 




05/07/26 




WFS90 


IRAS 


20126+4104 


20 


14 


25 


7 


+41 


13 


30 


04/05/14 


04/05/13 


05/07/21 


WFS91 


IRAS 


20205+3948 


20 


22 


20 





+39 


58 


21 


04/05/14 


04/05/13 


05/07/21 


WFS92 


IRAS 


20205+3948 


20 


22 


24 


9 


+39 


57 


55 




05/07/26 




WFS93 


IRAS 


20216+4107 


20 


23 


23 


9 


+41 


17 


42 




05/07/26 




WFS 94 


IRAS 


20293+3952 


20 


31 


12 


9 


+40 


03 


21 




05/07/27 




WFS95 


IRAS 


20319+3958 


20 


33 


49 


4 


+40 


08 


32 




04/05/15 




WFS96 


IRAS 


20332+4124 


20 


34 


58 


7 


+41 


34 


46 




04/05/15 




WFS97 


IRAS 


20332+4124 


20 


35 


01 


1 


+41 


34 


59 




04/05/15 




vv r oyo 


IRAS 




9(1 




uo 




4-4-1 




44 




f\K /n7 /97 




WFS99 


IRAS 


20343+4129 


20 


36 


06 


3 


+41 


39 


59 


04/05/17 


04/05/23 




WFS100 


IRAS 


20343+4129 


20 


36 


08 


1 


+41 


39 


58 


04/05/17 


04/05/23 




WFS101 


IRAS 


22134+5834 


22 


15 


08 


9 


+58 


49 


08 




05/07/23 




WFS102 


IRAS 


22551+6221 


22 


57 


04 


3 


+62 


37 


44 




05/07/23 




WFS103 


IRAS 


22551+6221 


22 


57 


07 


4 


+62 


37 


29 




05/07/23 




WFS104 


IRAS 


22551+6221 


22 


57 


11 


6 


+62 


36 


46 




05/07/23 




WFS107 


IRAS 


22570+5912 


22 


59 


05 





+59 


28 


23 


04/05/13 


04/05/13 


05/07/26 


WFS108 


IRAS 


23033+5951 


23 


05 


24 


8 


+60 


08 


14 


04/05/13 


04/05/13 


05/07/26 


WFS109 


IRAS 


23139+5939 


23 


16 


09 


8 


+59 


55 


31 


04/05/13 


04/05/13 


05/07/26 


WFS110 


IRAS 


23151+5912 


23 


17 


20 


4 


+59 


28 


51 


04/05/14 


04/05/22 


05/07/26 


WFS111 


IRAS 


23545+6508 


23 


57 


02 


1 


+65 


24 


38 


04/05/14 


04/05/22 


05/07/26 


WFS112 


IRAS 


23545+6508 


23 


57 


06 


4 


+65 


24 


49 


04/05/14 


04/05/22 


05/07/26 



Table 8. The telescope pointing coordinates are given along with the a s sociat ed IRAS source found within the frame. 
The WFS labels refer to the submm peaks identified in lWilliams et ah ( 20041 ). 



